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Abstract 

Searches for a scalar top quark and a scalar bottom quark have been performed using 
a total data sample of 56.8 pb _1 at a centre-of-mass energy of ps =183 GeV collected 
with the OPAL detector at LEP. No candidate events were observed. Combining this 
result with those obtained at lower centre-of-mass energies, the 95% C.L. lower limit on 
the scalar top quark mass is 85.0 GeV if the mixing angle between the supersymmetric 
partners of the left- and right-handed states of the top quark is zero. The lower limit is 
81.3 GeV, even if the scalar top quark decouples from the Z° boson. These limits were 
obtained assuming that the scalar top quark decays into a charm quark and the lightest 
neutralino, and that the mass difference between the scalar top quark and the lightest 
neutralino is larger than 10 GeV. The complementary decay mode of the scalar top quark 
in which it decays into a bottom quark, a charged lepton and a scalar neutrino was also 
studied. From a similar analysis, a mass limit on the light scalar bottom quark was 
set at 82.7 GeV for a mass difference between the scalar bottom quark and the lightest 
neutralino larger than 7 GeV, and at 84.0 GeV for the mass difference larger 10 GeV and 
the lightest neutralino heavier than 30 GeV. These limits were obtained assuming that 
the scalar bottom quark decays into a bottom quark and the lightest neutralino, and that 
a mixing angle between the supersymmetric partners of the left- and right-handed states 
of the bottom quark is zero. 
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1 Introduction 

Supersymmetric (SUSY) extensions JIJ of the Standard Model predict the existence of the 
bosonic partners of all known fermions. The scalar top quark (t), which is the bosonic partner 
of the top quark, can be the lightest charged supersymmetric particle for the following two 
reasons |2|. Firstly, one-loop radiative corrections to the t mass through Higgsino-quark loops 
and Higgs-squark loops are always negative. The correction is large for a heavy top quark mass 
as measured by the CDF and DO Collaborations [J3]. Secondly, the supersymmetric partners 
of the right-handed and left-handed top quarks (t R and ty) mix, and the resultant two mass 
eigenstates (ty and t 2 ) have a mass splitting. This mass splitting is expected to be very large 
due to the large top quark mass. It is possible that the lighter mass eigenstate (ty) can be 
lighter than any other charged SUSY particle, and also lighter than the top quark itself. The 
ty is the mixed state of t R and ty, i.e. ty = ty cos 6~ x + t R sin#p where 6~ t is a mixing angle. All 
SUSY breaking parameters [|lj are hidden in the 6~ t and a mass of ty. 

The scalar bottom quark (b) can also be light if tan/3, the ratio of vacuum expectation 
values of the two Higgs doublets, is larger than approximately 40. In this case, analogous 
mixing between the supersymmetric partners of the right- and left-handed states of the bottom 
quark (b R and ty) becomes large, and the resultant two mass eigenstates (bi and b 2 ) also have 
a large mass splitting [|4j. The mass of the lighter mass eigenstate (bi) may therefore be within 
the reach of LEP2. 

Scalar top quark pairs and scalar bottom quark pairs are produced in e + e _ annihilation 
via a virtual Z° boson or a virtual photon. In this paper it is assumed that R-parity [||] is 
conserved and that either y? or h is the on ly SUSY particle which is lighter than ty (bi), where 
y? is the lightest neutralino and z> is the scalar neutrino. The dominant decay mode of the ty 
with the above assumptions is expected to be eitherf] ty —> ex? or ty —> b u£ + . Both of these 
decay modes have been searched for. The dominant decay mode of the bi is expected to be 
bi —> by?. Under the assumption of R-parity conservation, the y? and u are invisible in the 
detector. Thus, tyty and bibi events are characterised by two acoplanar jets[] or two acoplanar 
jets plus two leptons, with missing energy. 

1 Through out this paper, all references to particle or decay implicitly include charge conjugation. 

2 Two jets not back-to-back with each other in the plane perpendicular to the beam axis. 
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The DO Collaboration has reported a lower limit [i6j on the ty mass of about 85 GeV (95% 
C.L.) for the case of ty —> cxi and that the mass difference between ty and Xi is larger than 
about 35 GeV. Searches at e + e“ colliders are sensitive to a smaller mass difference and mass 
limits for the ty have been obtained around the Z° peak (LEP1) assuming ty —> cxi■ A 95% 
C.L. lower limit of about 45 GeV was obtained for a mass difference larger than 5 GeV [;7]. 
Previous searches at centre-of-mass energies of \fs= 130, 136 0, 161 0 and 171 GeV 
improved the limit on the mass of the % to 66.8 GeV if the Q~ t is smaller than 7t/4 and if the 
mass difference between % and Xi is larger than 10 GeV. 

In 1997 the LEP e + e~ collider at CERN was run at centre-of-mass energies of 181-184 GeV. 
The luminosity weighted mean centre-of-mass energy was 182.7 GeV. In this paper direct 
searches for ty and bi using the data collected with the OPAL detector at these centre-of-mass 
energies are reported. The results shown here have been obtained by combining the results 
obtained at these new centre-of-mass energies with those obtained at y/s = 130, 136, 161 and 
171 GeV @,0. 

The phenomenology of the production and decay of the ty (bi) is described in section 2 
of the previous publication []TU|. In this paper, the OPAL detector and the event simulation 
for signal and background processes are given in section 2. In section 3, the data analysis is 
described and the results are given in section 4. 


10 , 11 


2 The OPAL Detector and Event Simulation 


2.1 The OPAL Detector 


The OPAL detector, which is described in detail in ref. |12|, is a multipurpose apparatus 
having nearly complete solid angle coverage. The central detector consists of a silicon strip 
detector and tracking chambers, providing charged particle tracking for over 96% of the full 
solid angle, inside a uniform solenoidal magnetic field of 0.435 T. A lead-glass electromagnetic 
(EM) calorimeter located outside the magnet coil is hermetic in the polar anglcf] range of 
| cos^| < 0.82 for the barrel region and 0.81 < | cos 6*| < 0.984 for the endcap region. The 
magnet return yoke consisting of barrel and endcap sections along with pole tips is instrumented 
for hadron calorimetry (HCAL) in the region |cos0| < 0.99. Four layers of muon chambers 
cover the outside of the hadron calorimeter. Calorimeters close to the beam axis measure the 
luminosity using small angle Bhabha scattering events and complete the geometrical acceptance 
down to 24 mrad. These include the forward detectors (FD) which are lead-scintillator sandwich 
calorimeters and, at smaller angles, silicon tungsten calorimeters (SW) |L3[ located on both 
sides of the interaction point. Tungsten shields were installed around the beam pipe in front 
of the SW detectors to reduce the amount of synchrotron radiation seen by the detector. The 
presence of the shield results in a gap in the SW acceptance between the polar angles of 28 and 
31 mrad. The gap between the endcap EM calorimeter and the FD is filled by an additional 
lead-scintillator electromagnetic calorimeter, called the gamma-catcher. 

3 A right-handed coordinate system is adopted, where the rr-axis points to the centre of the LEP ring, and 
positive z is along the electron beam direction. The angles 9 and (f> are the polar and azimuthal angles, 
respectively. 
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2.2 Monte Carlo Event Simulation of ty and bi 


Monte Carlo simulation of the production and decay of the ti was performed as follows 114 . 


The tyty pairs were generated taking into account initial-state radiation [jnj. The hadronisation 
process was subsequently performed to produce colourless ty-hadrons and other fragmentation 
products according to the Lund string fragmentation scheme (JETSET 7.4) 0, 0- The 
parameters for perturbative QCD and fragmentation processes were optimised using event 
shape distributions of the hadronic Z° decays measured by OPAL [17 . For the fragmentation 


of the ti, the fragmentation function proposed by Peterson et al. [0, 0 was used, where the 
parameter eq was set to 


eq = e b ml/mg (e b = 0.0038 [17], m h = 5 GeV) 


( 1 ) 


The ty-hadron was formed from a ty-quark and a spectator anti-quark or cliquark JT9 
The fragmentation products excluding the ty-hadrons carry less than 2% of the centre-of-mass 
energy. For the ty decaying into cxi, a colour string was stretched between the charm quark 
and the spectator. This colour singlet system was hadronised using the Lund scheme 0 


Gluon bremsstrahlung (QCD parton showering) was allowed in this process, and the Peterson 
function was also used for the charm quark fragmentation, where e c was set to 0.031 [0. The 
signals for the decays ty —> b £ + i> and bi —> by° were simulated in a similar manner. 


cXi, 


One thousand events were generated at each of 56 combinations of (my,m^o) for t b 
40 combinations of (rni^ma) for ty —> b £ + u and 40 combinations of for ty —> br + z>. 

The bibi events were generated for 48 combinations of (m bl ,m^o). The mixing angles of the 

ty and bi were set to zero when these events were generated. The dependence of the detection 
efficiencies on these mixing angles will be discussed in section 4.1. The generated events were 
processed through the full simulation of the OPAL detector 
chain was applied to the simulated events and the data. 


20 |, and the same event analysis 


2.3 Monte Carlo Event Simulation of Background Processes 

The background processes were simulated as follows: 

• Multijet hadronic events e + e _ —> qq(y) in which one or two jet momenta are mismeasured 

are a background for the large Am region (= m~ tl — rryo, my — m^, or rriy n — o). The 

PYTHIA [0 generator was used to simulate hadronic events. 

• r pairs, in which one of the r lepton decays into a low momentum electron and energetic 
neutrinos, are a background to acoplanar two-jet events. The KORALZ event generator [ |2I[ 
was used for the generation of r + r _ ( 7 ) and /r + /i - ( 7 ) events. The BHWIDE program [0 was 
used for e + e~ —> e + e _ ( 7 ) events. 

• Two-photon processes are the most important background for the case of a small mass differ¬ 
ence Am, since such signal events have small visible energy and small transverse momentum 
relative to the beam direction. Using the Monte Carlo generators PHOJET 0], PYTHIA [0 
and HERWIG [24], hadronic events from two-photon processes were simulated in which the 
invariant mass of the photon-photon system (M 77 ) was larger than 5.0 GeV. Monte Carlo 
samples for four-lepton events (e + e _ e + e _ , e + e _ /r + /r~ and e + e _ r + r _ ) were generated with the 


Vermaseren program 0. 
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• Finally, four-fermion processes in which at least one of the fermions is a neutrino constitute a 
serious background. The dominant contributions come from W + W - or y*Z 0 events. Since the 
interference effects of various diagrams are important, the grc4f generator |26j was used, which 
takes into account all interfering diagrams and includes initial-state photon radiation. 

These background events were also processed through the full detector simulation and the 
same event analysis chain as used for the data was applied. 


3 Analysis 

The present analysis is based on the data collected during the 1997 run of LEP2. Since the 
event topologies of ty —> cy'5 and bx —> by? are similar, the same selection criteria were used for 
both these modes (section 3.1 analysis A). In section 3.2 (analysis B), the selection criteria for 
ti —> b£ + z/ are discussed. These analyses are similar to those in ref. [jTp], and the same quality 
criteria as in ref. [|10j were used to select tracks and clusters. 

Variables used for the cuts, such as the total visible energy, E vis , the total transverse mo¬ 
mentum and the acoplanarity angle (defined below) were calculated as follows. First, the four- 
momenta of the tracks and those of the EM and HCAL clusters not associated with charged 
tracks were summed. Whenever a calorimeter cluster had associated charged tracks, the ex¬ 
pected energy deposited by the tracks was subtracted from the cluster energy to reduce double 
counting. If the energy of a cluster was smaller than the expected energy deposited by the 
associated tracks, the cluster energy was not used. Hadron calorimeter clusters were also used 
in calculating event variables. A large momentum-unbalance is occasionally caused by the fluc¬ 
tuation in the energy measurement of clusters in the hadron calorimeter because of the limited 
energy resolution. Therefore the transverse momentum and the visible mass calculated without 
the HCAL clusters were also used to reduce two-photon background processes. 

The following preselection criteria (PI - P2), which are common to both analyses A and B, 
were applied first. The numbers of events remaining after each cut are listed in Table |lj. For 
comparison, the table also shows the corresponding numbers of simulated events for background 
processes and for three samples of the simulated titi (ty —> cy?) and bibi events. 

(pi) The number of charged tracks was required to be at least four. The ratio of the number 
of good tracks to the total number of reconstructed tracks was required to be greater than 
0.2 to reduce beam-gas and beam-wall background events. The visible mass of the event, 
excluding the hadron calorimeter, was also required to be larger than 3 GeV. 

(P2) To reduce the background from two-photon processes and multihadronic events, where 
a jet axis was close to the beam direction, the total energy deposited had to be less than 
5 GeV in each SW detector, less than 2 GeV in each FD detector and less than 5 GeV in 
each side of the gamma-catcher. 

3.1 Analysis A: ty —» and bi —> 

The experimental signature for tyty(ty —> cy) 1 ) events and bibi events is an acoplanar two-jet 
topology with a large transverse momentum with respect to the beam axis. The fragmentation 
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functions of G and bi are very hard and the invariant mass of charm (bottom) quark and the 
spectator quark is small, therefore the jets are expected to be narrow. 

The event selection criteria are described below. The numbers of events remaining after 
each cut are also listed in Table [Tj. 

(Al) The visible energy in the region of |cosA| > 0.9 was required to be less than 20% of 
the total visible energy. In addition, the polar angle of the missing momentum direction, 
was required to satisfy |cosA miss | < 0.9 to reduce the two-photon and beam-gas 
events. 

(A2) Events from two-photon processes were largely removed by demanding that the event 
transverse momentum excluding the hadron calorimeter, P t , be greater than 4.5 GeV and 
that the transverse momentum including the hadron calorimeter, P f HCAL ; be greater than 
4.5 GeV. Fig. 1 shows the distribution of P t just before these cuts. 

(A3) The number of reconstructed jets was required to be exactly two. Jets were recon¬ 
structed using the Durham algorithm |27j with the jet resolution parameter of y cut 
= 0.005(E' vis / v / i) _1 . This FUg-dependent y cut parameter was necessary for good jet- 
reconstruction over a wide range of mr tl , rriy n and o. Fig. 2 shows the number of 
reconstructed jets before this cut. Furthermore, both reconstructed jets were required to 
contain at least two charged particles to reduce the t + t~ background where at least one 
of the r’s decayed into only one charged particle. 

(A4) The acoplanarity angle, 0 acop , is defined as 7r minus the azimuthal opening angle between 
the directions of the two reconstructed jets. To ensure the reliability of the calculation 
°f (ft&cop, both jet axes were required to satisfy | cos Ajetl < 0.95, where Aj et is the polar 
angle of the jet. The value of </> acop was required to be greater than 20°. Fig. 3 shows the 
distributions of </> acop just before this selection. 

(AS) ‘Softness’ was defined as where Mi and M 2 are the invariant masses of the 

two reconstructed jets, and E\ and E 2 are the energies of the jets. The signal events have 
low values of ‘Softness’, on the other hand, the two-photon events, which pass through 
the acoplanarity cut, have relatively large values. It was required that 1.5 x Softness < 
(P t — 4.5), where P t is given in units of GeV. Fig. 4 shows the scatter plots of Softness 
versus P t for data, the simulated two-photon events, the tyty and bibi events. 

(A6) The arithmetic mean of the invariant masses of the jets, Mj et , was required to be smaller 
than 8 GeV. When the invariant mass of the event, M v i S , was larger than 65 GeV, a 
harder cut M- ]et < 5 GeV was applied to reduce background from Wez/ events. Fig. 5 
shows the scatter plots of M vis versus Mj et for data, the simulated four-fermion events, 
the tit 1 and bjbi events. As shown in Fig. 5(c) and (d), jets from the ty and bi are 
expected to have low invariant masses, because the fragmentation function of the ty is 
hard and only a few particles are emitted from the fragmentation process of tit 1 . 

After all the cuts, no events were observed in the data, which is consistent with the expected 
number of background events of 2.0. The four-fermion processes are the dominant background 
processes. Uncertainties in the number of background for these processes will be discussed in 
Section 4.2. 
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data 

total 

bkg. 

qq(y) 

£ + £~(x) 

‘77’ 

4-f 

titj 

and bjbi 

m- tl (GeV) 







80 

80 

- 

m 6l (GeV) 







- 

- 

80 

iriy .o (GeV) 







75 

60 

60 

Presel. 1 

275501 

250435 

5028 

1229 

243137 

1042 

874 

967 

994 

Presel. 2 

141619 

124394 

3618 

1165 

118753 

859 

848 

933 

960 

cut (Al) 

22927 

19583 

1126 

236 

17651 

569 

731 

768 

809 

cut (A2) 

1441 

1461 

771 

156 

47.1 

488 

458 

745 

778 

cut (A3) 

362 

359 

298 

10.4 

16.3 

34.0 

305 

649 

717 

cut (A4) 

15 

25.3 

0.09 

0.06 

3.41 

21.8 

289 

582 

638 

cut (A5) 

15 

22.3 

0.08 

0.06 

0.38 

21.7 

229 

582 

636 

cut (A6) 

0 

1.97 

0.04 

0.04 

0.27 

1.63 

229 

576 

609 


Table 1: The remaining numbers of events normalised to the integrated luminosity of the data 
for various background processes are compared with data after each cut. Numbers for three 
simulated event samples of titi and bibi are also given (each starting from 1000 events). Before 
the cut (A2) was applied, there is a discrepancy between data and the simulated background 
processes, since the numbers of events expected from two-photon processes do not include the 
region M 77 < 5 GeV. 


The efficiencies for titi and bibi events are listed in Tables and |3]. Both efficiencies are 
30-60%, if the mass difference between the ty(bi) and Xi is larger than 10 GeV. A modest 
efficiency of about 20% is also obtained for a mass difference of 5 GeV for tyti events. In 
addition to effects included in the detector simulation, an additional efficiency loss of 3.6% 
(relative) arose from beam-related background in the SW, FD and gamma-catcher detectors 
estimated using random beam crossing events. The efficiencies given in Tables 1-5 do not 
include this correction, but it is included when deriving the mass limits. The efficiency at an 
arbitrary point was estimated using a polynomial fit to the efficiencies determined using the 
Monte Carlo simulations. 


3.2 Analysis B: ti —* b tv 

The experimental signature for titi (ti —> b tv) events is an two acoplanar jets plus two leptons 
with missing transverse momentum with respect to the beam axis. The momenta of leptons and 
the missing transverse momentum depend strongly on the mass difference between ty and v. To 
obtain optimal performance, two sets of selection criteria (analyses B-L and B-H) were applied 
depending on this mass difference. If the mass difference is smaller than or equal to 10 GeV, 
the momenta of leptons and the missing transverse momentum are relatively small. In such 
cases, it is difficult to identify leptons effectively and the dominant background comes from two- 
photon processes. Selections were optimised to reduce these two-photon events. When the mass 
difference is larger than 10 GeV, the momenta of leptons and the missing transverse momentum 
become large. In such cases the four-fermion processes become the dominant background 
processes. 

The numbers of events remaining after each cut are listed in Tables [| and [5]. For comparison, 

























m- u (GeV) 

50 

55 

60 

65 

70 

75 

80 

85 

Am 









3 GeV 

11 

9 

7 

6 

5 

3 

2 

2 

5 GeV 

29 

30 

27 

28 

27 

25 

23 

22 

10 GeV 

40 

44 

45 

48 

49 

48 

48 

49 

20 GeV 

36 

45 

49 

52 

53 

58 

58 

60 

mj 2 

35 

38 

41 

45 

44 

47 

51 

52 

mq - 10 GeV 

27 

28 

28 

29 

32 

35 

40 

40 


24 

27 

29 

29 

33 

33 

37 

36 


Tabic 2: The detection efficiencies in percent for titq, in which ty decays into cy? for different ty 
masses and Am values, where Am is mq — m^o . The statistical fluctuations of these efficiencies 
are about 2% (absolute). 


(GeV) 

50 

55 

60 

65 

70 

75 

80 

85 

Am 









7 GeV 

34 

34 

36 

36 

35 

38 

37 

36 

10 GeV 

40 

43 

46 

46 

50 

49 

50 

49 

20 GeV 

36 

41 

48 

51 

55 

59 

61 

63 

"ffii/ 2 

29 

34 

35 

37 

37 

41 

47 

50 

mg 1 - 10 GeV 

22 

22 

22 

25 

25 

27 

32 

33 

m bi 

21 

23 

21 

22 

23 

26 

31 

31 


Tabic 3: The detection efficiencies in percent for bibi for different bi masses and Am val¬ 
ues, where Am = m^ i — rri^o. The statistical fluctuations of these efficiencies are about 2% 
(absolute). 

the tabic also shows the corresponding numbers for simulated events for background processes 
and for two samples of simulated titi events, in which the branching fraction to each lepton 
flavour is assumed to be the same. 


3.2.1 Small mass difference case 


For the case of small mass difference (Am < 10 GeV), the following selections were applied. 


(B-Ll) The visible energy in the region of | cosd| > 0.9 was required to be less than 10% of 
the total visible energy. In addition, | cos fl mi SS | < 0.8 was required. 

(B-L2) Both P t and p t HCAL were required to be greater than 5 GeV. 


(B-L3) The number of charged tracks was required to be at least six. Furthermore, the 
number of reconstructed jets was required to be at least four, because the signal should 
contain two hadronic jets plus two isolated leptons. Jets were reconstructed using the 
Durham algorithm [^7J with the jet resolution parameter of y cut = 0.004. Fig. 6 shows 
the distributions of the number of reconstructed jets before this selection. 


(B-L4) To examine the acoplanarity of the events, jets were reconstructed using the Durham 
algorithm where the number of jets was forced to be two. To ensure a good measurement 


9 









































of acoplanarity angle, | cos 0jet| < 0.95 was required for both reconstructed jets. Finally, 
the acoplanarity angle between these two jets was required to be greater than 15°. In the 
three-body decay, the transverse momentum carried by the z> with respect to the original 
ti-momentum is smaller than that of x± i n the two-body decay. When the tq is light, the 
outgoing z> is strongly boosted toward the direction of the parent tq. Hence 0 acop for the 
signal becomes small. This is the reason for the use of a looser acoplanarity angle cut. 

(B-L5) The total visible energy normalised to the centre-of-mass energy, E v j S / a /s, was required 
to be smaller than 0.2 to reject four-fermion events. 


No events were observed in the data after the above cuts. The number of expected back¬ 
ground events is 1.1. Two-photon processes are the dominant background. Uncertainties in the 
number of expected background events will be discussed in Section 4.2. 



data 

total 

bkg. 

qq(y) 

£+r(y) 

‘77’ 

4-f 

titi 

m- tl (GeV) 







80 

80 

mo (GeV) 







73 

70 

cut (B-Ll) 

11713 

9354 

904 

134 

7859 

457 

717 

695 

cut (B-L2) 

1064 

1076 

592 

78.3 

18.2 

388 

139 

440 

cut (B-L3) 

287 

295 

94.8 

0.02 

1.84 

198 

137 

435 

cut (B-L4) 

60 

62.6 

3.33 

0.00 

0.97 

58.3 

118 

379 

cut (B-L5) 

0 

1.07 

0.01 

0.00 

0.97 

0.09 

118 

378 


Table 4: The remaining numbers of events normalised to the integrated luminosity of the 
data for various background processes are compared with data after each cut. Numbers for 
two simulated event samples of tqti are also given (each starting from 1000 events). In these 
samples, the branching fraction to each lepton flavour is assumed to be the same. Before the cut 
(B-L2), there is a discrepancy between data and the simulated background processes, since the 
numbers of events expected from two-photon processes do not include the region M 77 < 5 GeV. 


3.2.2 Large mass difference case 

The selections for a large mass difference (Am > 10 GeV), are described below. Because the 
momenta of the leptons for this case are high enough to be identified, it was required that events 
contained at least one lepton. Then the other cuts (B-Hl and B-H3) were relaxed compared to 
the small mass difference case. 

(B-Hl) Cut (Al) was applied to reduce two-photon and beam-gas events. 

(B-H2) Cut (B-L2) was applied to reject two-photon events. 

(B-H3) The number of charged tracks was required to be at least six, and the number of 
reconstructed jets was required to be at least three. Jets were reconstructed with the jet 
resolution parameter of y cut = 0.004. 

(B-H4) Cut (B-L4) was applied to reject multihadronic events. 
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(B-H5) A candidate event was required to contain at least one lepton. Leptons were identified 
in the following way: electrons were selected if they satisfied either of the two identification 
methods described in Ref. [^8|], and muons were identified using the two methods described 
in ref. [^9|. The track momentum of the electron or muon candidate was required to be 
larger than 2 GeV. A jet reconstructed in cut (B-H3) was identified as a tau decay if it 
contained only one or three charged tracks, the invariant mass of the charged particles 
in the jet was smaller than 1.5 GeV, the invariant mass including energies deposited in 
the calorimeters was smaller than 2 GeV and the scalar sum of momenta of the charged 
tracks was larger than 2 GeV. 


(B-H6) The invariant mass excluding the most energetic lepton was required to be smaller 
than 60 GeV in order to reject W + W” —► i/£qq' events. As shown in Fig. 7, a large 
fraction of four-fermion events was rejected using this requirement. 


(B-H7) The visible mass of the events, M vis , must be smaller than 80 GeV to reduce W + W 
background events in which one of W ±! s decayed into tv and the other into qq (g). If one 
jet from qq (g) was misidentified as tau lepton, this event could pass through the previous 
cut (B-H6). Such events were rejected by this requirement. Fig. 8 shows the distribution 


No events were observed in the data after the above cuts. The number of expected back¬ 
ground events was 2.1. Uncertainties of the expected background events will be discussed in 
section 4.2. The detection efficiencies for titi events are listed in Table || The efficiencies of 
both selection criteria (B-L and B-H) are presented in this table. As shown in this table, the 
detection efficiencies for tq —> b r + v T are slightly smaller than the case in which the branching 
fraction to each lepton flavour is assumed to be the same. 



data 

total 

bkg. 

qq(y) 

FlrO) 

‘77’ 

4-f 

titi 

m~ tl (GeV) 







80 

80 

80 

mo (GeV) 







70 

60 

40 

cut (B-Hl) 

22927 

19603 

1127 

236 

17670 

567 

856 

792 

784 

cut (B-H2) 

1360 

1381 

724 

151 

32.2 

474 

490 

748 

776 

cut (B-H3) 

857 

915 

466 

1.48 

8.66 

438 

490 

743 

764 

cut (B-H4) 

209 

213 

15.3 

0.36 

2.60 

195 

424 

665 

660 

cut (B-H5) 

162 

162 

6.59 

0.33 

0.32 

155 

369 

639 

637 

cut (B-H6) 

4 

6.35 

0.17 

0.14 

0.22 

5.83 

369 

639 

606 

cut (B-H7) 

0 

2.05 

0.15 

0.06 

0.22 

1.62 

369 

639 

581 


Table 5: The remaining numbers of events normalised to the integrated luminosity of the 
data for various background processes are compared with data after each cut. Numbers for 
three simulated event samples of tity are also given (each starting from 1000 events). In these 
samples, the branching fraction to each lepton flavour is assumed to be the same. Before the cut 
(B-H2), there is a discrepancy between data and the simulated background processes, since the 
numbers of events expected from two-photon processes do not include the region M 77 < 5 GeV. 
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ti —■> b iv 

the equal branching fractions for £= e, /i, r 

m ti (GeV) 

50 

55 

60 

65 

70 

75 

80 

85 

Am (B-L) 









7 GeV 

9 

11 

12 

12 

12 

12 

12 

9 

10 GeV 

16 

22 

30 

32 

35 

39 

38 

37 

Am (B-H) 









10 GeV 

30 

33 

36 

37 

36 

40 

37 

35 

20 GeV 

45 

52 

56 

60 

61 

63 

64 

65 

mj 2 

41 

46 

48 

52 

50 

51 

58 

60 

mq-10 GeV 

25 

26 

27 

27 

29 

26 

28 

27 

ti —> brz> T , 100% branching fraction 

m- tl (GeV) 

50 

55 

60 

65 

70 

75 

80 

85 

Am (B-L) 









7 GeV 

6 

6 

8 

8 

7 

6 

6 

4 

10 GeV 

16 

22 

26 

29 

31 

32 

29 

28 

Am (B-H) 









10 GeV 

19 

20 

21 

21 

22 

23 

22 

19 

20 GeV 

36 

40 

42 

46 

49 

50 

51 

52 

™%J 2 

35 

38 

44 

44 

48 

50 

54 

51 

m - tl -10 GeV 

26 

29 

31 

29 

33 

32 

32 

35 


Table 6: The detection efficiencies in percent for titi, in which tq decays into b£z> (£ = e, n, r). 
The upper half of the table shows the case in which the branching fraction to each lepton 
flavour is the same and the lower half shows the worst case in which the branching fraction of 
ti —»■ b tu t is 100%. In both tables, Am is defined as m~ tl — m„. The efficiencies in the first two 
lines of each table were obtained using the analysis B-L, the the last four lines using analysis 
B-H. 


4 Results 

No evidence for tqti and bibi pair-production has been observed in the data. The data are 
consistent with the expected background of 1.9 events in analysis A, and 1.0 and 2.0 events^] 
in analysis B for an integrated luminosity of 56.8 pb _1 . The sum of the expected number of 
background is 4.5 events subtracting overlap between three analyses. Uncertainties of expected 
background will be discussed in section 4.2. 

4.1 Systematic Errors in the number of expected signal events 

The following sources of systematic error on the expected number of the signal events were 
taken into account: 


4 Tliese numbers were corrected for the inefficiency due to beam-related background events. 
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1. The statistical error of the signal Monte Carlo simulation is 2-10% depending on detec¬ 
tion efficiencies. 


2. The dependence of the detection efficiency on the mixing angle: 

The energy distribution of the initial-state radiation depends on the mixing angle of the 
ty (bi), because it influences the coupling between the ty (bi) and the Z°. When the 
coupling is large, the initial-state radiation is hard. The detection efficiencies therefore 
depend on 6~ t (0g). However, the detection efficiencies in Tables |j. |3| and || were calculated 
using the simulated events which were generated for 6~ t = 9^ = 0.0. 

The detection efficiencies in the two extreme cases of ty decoupled from the Z° (0i = 0.98) 
and ty = ty (&i = 0.0) were compared for various my values. The difference was always 
found to be within 2%. The effect on efficiencies for ty —> biv and bi —> by? was also 
checked and similar results were obtained. The systematic error due to the dependence 
on the mixing angle was taken to be 2%. 


3. Fragmentation function for ty: 

The fragmentation scheme proposed by Peterson et al. was used, with the fragmentation 
parameter ey determined by equation (1). The error on ey was propagated from = 

± 0.26 [[nj and dm^/m^ = ± 0.06 [[I0j, corresponding to <5ey/ey = ± 0.27. The systematic 
error in the efficiencies due to this uncertainty was evaluated by altering the ey parameter 
by one standard deviation for several combinations of (my, m^o) and (my, ms). For the 


ty —> cy® mode, the detection efficiencies changed by no more than 5% over the my range. 
The relative changes for the ty —> b£v mode were found to be 6-10%, and they depended 
mainly on my. 

To estimate the dependence on the fragmentation scheme, the fragmentation function 
proposed by Bowler |3l]| was used, because the shape of this fragmentation function is very 
different from that of the Peterson function. The relative difference in efficiencies between 
the two fragmentation models was 2-3% for the ty —> cyj mode, which was smaller than 
that due to the variation of the ey parameter used in the Peterson fragmentation scheme. 
The systematic error due to the dependence on the fragmentation model was taken to be 
3%. For the ty —> b£u mode, the relative difference in efficiencies was found to lie between 
4-8%, where the range was mainly due to my. 


4. Fragmentation function for by 

The error due to the fragmentation function for bi was also estimated using the meth¬ 
ods described above. The uncertainty in eg made a relative difference of 4-6% in the 
efficiencies. 


5. Fragmentation of the charm and bottom quarks: 

The error in the efficiencies for the iy —> cyj mode, due to the uncertainty in e c , was 


estimated to be typically 3% by changing e c by 5e c /e c = ±0.35 |L7 


The uncertainty in the 6b parameter also contributes to the error in the efficiencies for 
the ty —> b£u and bi —» by? modes. As mentioned above, the eb parameter was simul¬ 
taneously changed by ±26% when ey and eg x were altered. The systematic error due to 
the uncertainty on 6b is therefore taken into account in the errors 5ey and 
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6. Fermi motion of the spectator quark in tq (bi) -hadron decay: 

Due to the Fermi motion of the spectator quark the invariant mass of the hadronic decay 
products of a tq (bi) -hadron varies. For tq —> cx° { and tq —> by° modes this effect is not 
negligible when Am is large. The systematic error in the efficiencies due to the Fermi 
motion was evaluated by altering the mass of spectator quarks. For the case of a 80 GeV 
tq (bi) and a massless neutralino the efficiency varies by up to 6% (8%) due to the jet 
mass cut (A6). 

7. Lepton identification: 

The systematic error on electron identification was estimated to be 4% and the error 
on muon identification was 2%. The systematic error on tau identification is dominated 
by the uncertainties in the fragmentation of the bottom quark, which has already been 
included in the uncertainty in the 6b parameter. A conservative error of 4% was applied 
for all types of leptons. 

8. Systematic errors due to imperfections in the Monte Carlo simulation of P t , P f HCAL ; the 
number of reconstructed jets, E vis and M vis were estimated to be 3%. 

9. The integrated luminosity was calculated using the SW detector. The systematic error 
on this luminosity was determined to be 0.26% (stat.) and 0.41% (syst.). 

10. The systematic error due to the uncertainty on the trigger efficiency was estimated to be 
negligible. This is expected because of the requirement of at least four good tracks. 


The various systematic errors are summarised in Table [7j. These systematic errors were 
considered to be independent and the total systematic error was calculated as the quadratic 
sum of the individual errors. These systematic errors were treated as in Ref. 
the limits. 


32] in calculating 


Sources 

U -»■ cxi 

tq —■> b £i> 

bi -»• by? 

Statistical error of MG 

2-10% 

6- t dependence 

2% 

2% 

- 

6 £ dependence 

- 

- 

2-4% 

Uncertainty on eq 

5% 

6-10% 

- 

Uncertainty on 

- 

- 

4-6% 

Fragmentation scheme 

3% 

4-8% 

4-10% 

Uncertainty on e c 

3% 

- 

- 

Uncertainty on eb 


Included in the 
uncertainties of e~ tl and eg, 

Spectator Fermi motion 

3-6% 

4% 

3-8% 

Uncertainty of lepton ID 

- 

4% 

- 

Detector simulation 

3% 

Luminosity 

0.5% 

Trigger efficiencies 

negligibly small 


Table 7: The summary of the systematic errors on the expected number of the signal events. 
The range of these errors depend on the mass of tq and bi. 
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4.2 Systematic Errors in the number of background events 

The two-photon and four-fermion processes are the dominant background. Systematic errors 
in the expected number of these processes are discussed here. Since no event was observed in 
data, the expected number of background events was not subtracted to calculate limits. Then 
these errors were not used to calculate limits. 


4.2.1 Two-photon processes 

The systematic errors are mostly dominated by the Monte Carlo statistics for two-photon 
processes. The statistical fluctuations of the expected numbers (relative errors) are ±0.23(85%), 
±0.49(50%) and ±0.22(100%) for analysis A, B-L and B-H, respectively. 

Furthermore the uncertainty on the modelling of the two-photon processes was checked with 
data. In order to select two-photon events the visible energy was required to be smaller than 
20% of y/s, the charged multiplicity to be at least four, the visible invariant mass to be larger 
than 5 GeV and the forward detector vetoes (cut P2) were required. The P t distributions of the 
selected events from data were compared with Monte Carlo. The shapes of the distributions 
agree with each other, but there is an uncertainty of 30% in the normalisation. 


4.2.2 Four-fermion processes 


Uncertainties in the generators of the four-fermion processes were estimated by comparing 

and PYTHIA 


grc4f with the Excalibur [33 


15] generators. The background events predicted 


by these different generators are summarised in Table 8. The differences larger than statistical 
fluctuations were found especially in analysis A. The difference between grc4f and Excalibur 
comes mainly from the region of M < 40 GeV for e + e~ —> ei/qq' process. On the other 


hand, the difference between grc4f and PYTHIA comes from the region of M qq > 15 GeV for 


e 1 e 


7 


—> qquP process. The prediction on this process by PYTHIA is about 60% of 
grc4f. These differences were considered as the systematic errors in the four-fermion processes; 
1.63 ± 0.13 (stat.) tgjg (sys.) for analysis A, 0.09 ± 0.03 (stat.) 1 q!o 6 (sys.) f° r analysis B-L 
and 1.62 ±0.13 (stat.) lo!ii (sys.) f° r analysis B-H. 



grc4f 

Excalibur 

PYTHIA 

analysis A 

1.63 ±0.13 

2.15 ±0.15 

1.30 ±0.10 

analysis B-L 

0.09 ±0.03 

0.21 ±0.05 

0.03 ±0.01 

analysis B-H 

1.62 ±0.13 

1.85 ±0.13 

1.44 ±0.13 


Table 8: The expected numbers of four-fermion background processes predicted by different 
three generators. The errors shown in this table are the statistical fluctuations. 

Events with two jets plus missing transverse momentum were checked with data to study 
the differences in analysis A. To select events with two jets and large P t coming from the 
four-fermion processes, the cuts PI, P2 and Al were applied, and the P t was required to be 
greater than 10 GeV to reject two-photon processes completely. Furthermore d\ 3 {= 2 / 23 -Evis 2 ) 
was required to be smaller than 50 GeV 2 to select clear two-jet events, where y 2 3 is the jet 
resolution parameter from 2 jets to 3 jets using the Durham algorithm. Finally, the acoplanarity 
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DATA 

grc4f 

Excalibur 

PYTHIA 

0 a cop > 10° 

19 

17.2 ±0.4 

18.9 ±0.4 

17.7 ±0.1 

0 acop > 100° 

2 

1.87 ±0.15 

3.25 ±0.18 

1.22 ±0.09 


Tabic 9: The remaining numbers of events in data and three different generators. The errors 
shown in this tabic are the statistical errors. 

angle, 0 acop , was required to be larger than 10° or 100°. After these selections, the observed 
numbers in data were compared with the predictions by these three generators (Table 9). 

In the region of </> acop > 100°, large differences were observed in the predictions of these 
three generators. But all three predictions are consistent with the data, since the statistics of 
data is too small. We need data with higher statistics to study the four-fermion generators. 


4.3 Mass Limits 

4.3.1 Scalar top quark ty 


To calculate mass limits, the number of signal events passing through the event selections is 
determined as a function of my, rn^o (or me) and 9~ x . 


Figs. 9(a), 10(a) and 11(a) show the 95% C.L. excluded regions in the ( 9 X , m~ tl ) plane for 
the ty —► cxi, ty —> b£z/ (£= e,p,r) and ty —> b tv decay modes, respectively. The branching 
fraction to each lepton flavour £ + depends on the composition of the lightest chargino [[10] . As 
the chargino becomes Higgsino-like, the branching fraction into b t + v t becomes large. In the 
limit that the chargino is the pure Wino state, the branching fraction to each lepton flavour 
is the same. Two extreme cases in which the branching fraction to each lepton flavour is the 
same, or the branching fraction into b t + v t is 100%, were considered here. The 95% C.L. mass 
bounds are listed in Table [H] for various values of 9 $. Assuming that the ty decays into cy5, an d 
that the mass difference between the ty and the is greater than 10 GeV, the ty is found to be 
heavier than 85.0 GeV, if 9 X = 0. A lower limit of 81.3 GeV is obtained even if the ty decouples 
from the Z° boson. When the ty decays into b£u, the lower limit on m Xl is 83.6 GeV, assuming 
that the mass difference between ty and v is greater than 10 GeV, that 9~ t = 0 and that the 
branching fraction to each lepton flavour is the same. The 95% C.L. excluded regions in the 
and {m~ tl , me) planes are shown in Figs. 9(b), 10(b) and 11(b) for various values of 


m 


'ti 


m xV 


9- t . In these figures, the limits calculated by the expected background events are superimposed 
to show the sensitivity of these analyses. 


4.3.2 Scalar bottom quark b x 

To calculate mass limits, the number of signal events passing through the event selections is 
determined as a function of ra^, m^o and 9 g. Fig. 12(a) shows the 95% C.L. excluded regions in 
the (#g, my) plane for the mass difference of A m(= my — m^ o) > 7 GeV and Am > 10 GeV. 
The numerical mass bounds are listed in Table [TT] for various 9 y The lower limit on the bi-mass 
is 82.7 GeV, if Am is greater than 7 GeV and 9 ^ = 0. The 95% C.L. excluded regions in the 
(my, m^o) plane are shown in Fig. 12(b) for various 9^. In this figure, the limit calculated by 
the expected background events is also superimposed to show the sensitivity of this analysis. 
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Lower limit on mq (GeV) 


% -*• ex ? 

ty — ► b £ z > 
l — e, n, t 

ti — * b tv t 

Br = 100% 

6~ t (rad) 

Am > 5 GeV 

Am > 10 GeV 

Am >10 GeV 

Am >10 GeV 

0.0 

81.2 

85.0 

83.6 

80.0 

< ±vr 

80.0 

84.2 

82.5 

79.0 

< 4 ^ 

76.8 

82.0 

79.7 

76.1 

0.98 

75.8 

81.3 

79.2 

75.0 


Table 10: The excluded my region at 95% C.L. (Am = m~ tl — rn^o or mq — mo). 


Lower limit on mg i (GeV) (bi —> by?) 

6 b ( rad ) 

Am >7 GeV 

Am >10 GeV 
mp > 30 GeV 

0.0 

82.7 

84.0 

<L 

81.0 

82.6 

<b 

71.9 

76.2 

1.17 

54.4 

63.7 


Table 11: The excluded m^ region at 95% C.L. (Am = m ]n — m^o) 

5 Summary and Conclusion 


A data sample collected using the OPAL detector corresponding to an integrated luminosity 
of 56.8 pb -1 at -y/s =183 GeV has been analysed to search for pair production of the scalar 
top quark and the scalar bottom quark predicted by supersymmetric theories. R-parity was 
assumed to be conserved. No events remained after the selection cuts. 

The 95% C.L. lower limits on the scalar top quark mass are 85.0 and 82.0 GeV, if the mixing 
angle of the scalar top quark is 0 and smaller than respectively. Even if the ty decouples 
from the Z° boson, a lower limit of 81.3 GeV is obtained. These limits were obtained assuming 
that the scalar top quark decays into a charm quark and the lightest neutralino and that the 
mass difference between the scalar top and the lightest neutralino is larger than 10 GeV. 

Assuming a relatively light scalar neutrino (37.1 GeV [[L|, [3DJ < rrip < m~ tl — mb) the 
complementary decay mode of the scalar top quark in which it decays into a bottom quark, a 
charged lepton and the scalar neutrino has also been studied. If the mass difference between 
the scalar top quark and the scalar neutrino is greater than 10 GeV and if the mixing angle of 
the scalar top quark is 0, the 95% C.L. lower limit on the scalar top quark mass is 83.6 GeV. 
This limit is obtained assuming that the branching fraction to each lepton flavour is the same. 

A mass limit on the light scalar bottom quark is found to be 82.7 GeV (95% C.L.), assuming 
that the mass difference between the scalar bottom quark and the lightest neutralino is greater 
than 7 GeV and that the mixing angle of the scalar bottom quark is zero. If the mass difference 
is greater than 10 GeV and the lightest neutralino is heavier than 30 GeV, a mass limit on the 
light scalar bottom quark is 84.0 GeV (95% C.L.). 
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Figure 1: The distributions of P t after cut (Al) for background (histograms) and data in 
(a), and for titi predictions in (b). The arrows in these figures show the selection criteria. 
In (a) the distribution of the data is shown by the points with error bars. The predictions 
from background processes are also shown: dilepton events (cross-hatched area), two-photon 
processes (grey area), four-fermion processes (singly-hatched area), and multihadronic events 
(open area), (b) shows predictions for tyti in which ty decays into ex?- The continuous line 
histogram is for (my, m^o) =(80 GeV, 60 GeV), and the dashed line is for (80 GeV, 75 GeV), 
starting from 1000 generated events each. 


OPAL 



(histograms) and data in (a), and for titi predictions in (b). The conventions for the various 
histograms are the same as in Fig. 1. 
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Figure 3: The distributions of the acoplanarity angle after cut (A3) for background (his¬ 
tograms) and data in (a), and for the titi predictions in (b). The conventions for the various 
histograms are the same as in Fig. 1. 


22 






























OPAL 



P t (GeV) P t (GeV) 


Figure 4: Scatter plots of ‘Softness’(see text) versus P t after cut (A4) for (a) data, (b) simulated 
two-photon processes, (c) the Monte Carlo simulation of titi signals for (m^m^o) = (80 GeV, 

75 GeV), and (d) the Monte Carlo simulation of bibi signals for (80 GeV, 73 GeV). 

The scale of P t in (a) is different from the other figures. The simulated events are not normalised 
to the luminosity. For the two-photon processes, the corresponding luminosity is 253 pb _1 . The 
event samples of titi and bibi each start from 1000 events. 
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Figure 5: Scatter plots of M vis versus M- ]et after cut (A5) for (a) data, (b) simulated four- 
fermion processes, (c) the Monte Carlo simulation of titi signals for (777^ ,771^.0) = (80 GeV, 

60 GeV), and (d) the Monte Carlo simulation of bjbi signals for ( 777 ^, 777 ^ 0 ) = (80 GeV, 60 GeV). 
The simulated events are not normalised to the luminosity. For the four-fermion processes, the 
corresponding luminosity is 5000 pb^ 1 . The event samples of tit 1 and bibi each start from 
1000 events. 


24 



OPAL 




in these figures show the selection criteria, (a) shows the distribution of the data with error 
bars. The predictions from background processes are also shown: dilepton events (cross-hatched 
area), two-photon processes (grey area), four-fermion processes (singly-hatched area), and the 
multihadronic events (open area), (b) shows predictions for titi in which ty decays into b£i). 
The continuous line histogram is for (m^, mo) =(80 GeV, 70 GeV), and the dotted line is for 
(80 GeV, 73 GeV). In these samples, the branching fraction to each lepton flavour is assumed 
to be the same. 
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Figure 7: 

(B-H5). The arrows in these figures show the selection criteria. The conventions for the various 
histograms in (a) are the same as in Fig. 6. In (b) the continuous line histogram is for (m^, 
mo) =(80 GeV, 60 GeV), and the dotted line is for (80 GeV, 40 GeV). In these samples, the 
branching fraction to each lepton flavour is assumed to be the same. 
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Figure 8: The distributions of the visible mass after cut (B-H6) for background (histograms) 
and data in (a), and for the tyti signal predictions in (b). The conventions for the various 
histograms are the same as in Fig. 7. 

OPAL (t, -> c x») 




Figure 9: The 95% C.L. excluded regions assuming that the ty decays into cxi■ (a) The 
excluded regions in the (9i, rn ^) plane for a mass difference Am (= — rn^o) > 10 GeV, and 

Am > 5 GeV. The cross-hatched region has already been excluded by the search at LEP1 [7|J. 
(b) The excluded regions in the (m^, m^o) plane, for a mixing angle of ty of 0.0 and 0.98 
rad. The solid lines show the actual limits, and the thin lines show the limits calculated only 
with the expected number of background events. The cross-hatched region has already been 
excluded by the search at LEP1 [0. The singly-hatched region has been excluded by the DO 
Collaboration [6]. The dash-dotted straight line shows the kinematic limit for the ty —> cy) 1 
decay. 
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Figure 10: The 95% C.L. excluded regions assuming that the ti decays into b£D and that 
the branching fraction to each lepton flavour is the same, (a) The excluded regions in the 
(6*y my) plane where the mass difference between the ty and the z> is greater than 8 or 10 GeV. 
The dash-dotted straight line shows the kinematic limit for this decay, since a v lighter than 

. (b) The excluded regions in the (my, m„) plane, for a 



34, 


37.1 GeV has been excluded 
mixing angle of the ty assumed to be 0.0 and 0.98 rad. The solid lines show the actual limits, 
and the thin lines show the limits calculated only with the expected number of background 
events. The dash-dotted horizontal line shows the limit on rrip obtained at LEP1, and the 
dash-dotted diagonal line shows the kinematic limit for the ty —> b£v decay. 
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Figure 11: The 95% C.L. excluded regions assuming that ty always decays into brz> T . (a) The 
excluded regions in the (0y my) plane where the mass difference between the ty and the v T 
is greater than 8 or 10 GeV. The dash-dotted straight line shows the kinematic limit for this 
decay, (b) The excluded regions in the (my, mp) plane, for a mixing angle of the ty assumed to 
be 0.0 and 0.98 rad. The solid lines show the actual limits, and the thin lines show the limits 
calculated only with the expected number of background events. The dash-dotted horizontal 
line shows the limit on rrip obtained at LEP1 |34| and the dash-dotted diagonal line shows 
the kinematic limit for the ty —> b tv decay. 
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Figure 12: The 95% C.L. excluded regions assuming that the bi decays into by?- 

(a) The excluded region in the (6% rn ^ i ) plane for a mass difference, Am (= — m^o), Am > 

10 GeV and yj is heavier than 30 GeV. The excluded region for Am > 7 GeV is also shown. 

(b) The excluded regions in the (m^, m^o) plane, for a mixing angle of the bi assumed to be 
0.0 and 1.17 rad. The solid lines show the actual limits, and the thin lines show the limits 
calculated only with the expected number of background events. 
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